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DNA methylation is evolutionarily conserved. Vertebrates exhibit high, widespread DNA 
methylation whereas invertebrate genomes are less methylated, predominantly within 
gene bodies. DNA methylation in invertebrates is associated with transcription level, 
alternative splicing, and genome evolution, but functional outcomes of DNA methylation 
remain poorly described in lophotrochozoans. Recent genome-wide approaches improve 
understanding in distant taxa such as molluscs, where the phylogenetic position, and life 
traits of Crassostrea gigas make this bivalve an ideal model to study the physiological 
and evolutionary implications of DNA methylation. We review the literature about DNA 
methylation in invertebrates and focus on DNA methylation features in the oyster. Indeed, 
though our MeDIP-seq results confirm predominant intragenic methylation, the profiles 
depend on the oyster's developmental and reproductive stage. We discuss the perspective 
that oyster DNA methylation could be biased toward the 5'-end of some genes, depending 
on physiological status, suggesting important functional outcomes of putative promoter 
methylation from cell differentiation during early development to sustained adaptation of 
the species to the environment. 
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DNA methylation is an epigenetic mechanism of great biologi- 
cal significance which is widely conserved in evolution (Zemach 
et al., 2010). DNA methylation in prokaryotes mostly affects 
adenines and is implicated in a broad range of processes from 
the control of DNA replication (review in Collier, 2009) to the 
defense against bacteriophages (Zaleski et al, 2005). Besides, bac- 
terial DNA methylation contributes to the inheritable control of 
gene expression, maintaining phenotypes through cell genera- 
tions, as observed for virulence genes (Heusipp et al., 2007). 

In contrast, DNA in eukaryote genomes is predominantly 
methylated on cytosines (Zemach et al., 2010). However, sig- 
nificant differences exist between kingdoms. In animal genomes 
cytosines are methylated mostly within CpG dinucleotides 
(Doskocil and Sorm, 1962; Lister et al., 2009), whereas plants 
(Cokus et al., 2008), and fungi (Antequera et al., 1985; Selker and 
Stevens, 1987) can also have methylcytosines in the CHG or CHH 
contexts. In addition, genomes display markedly different methy- 
lation profiles between vertebrates and invertebrates. Vertebrates 
exhibit a high DNA methylation, whereas invertebrate genomes 
are far less methylated (Suzuki et al, 2007; Suzuki and Bird, 2008; 
Feng et al, 2010; Zemach et al., 2010). 

In vertebrates, DNA methylation affects the whole genome, 
but eventually drops in the 5' regions of genes (Zemach et al., 
2010) which divide into two groups regarding their CpG content 
(Antequera, 2003). The low-CpG promoters are hypermethy- 
lated and drive the transcription of tissue-specific genes. Instead, 
the high-CpG promoters where CG dinucleotide-rich regions 



define "CpG islands" (CGIs) are hypomethylated and control 
widely expressed genes (Elango and Yi, 2008). Depending on 
their density, methylcytosines in promoters restrict the access 
of the transcription machinery to transcription start sites due 
to CpG-binding proteins recruitment (Antequera, 2003; Deaton 
et al., 2011; Vinson and Chatterjee, 2012) and/or the regula- 
tion of DNA/histone interactions within nucleosomes (Cedar and 
Bergman, 2009; Zentner and Henikoff, 2013) thereby diminishing 
transcription (Hsieh, 1994; Riviere et al, 2011). 

Invertebrates generally display a heterogenic or "mosaic" 
methylation profile with predominant methylation of transcrip- 
tion units (Suzuki et al., 2007; Elango et al., 2009; Gavery and 
Roberts, 2010; Walsh et al, 2010; Gadau et al, 2012; Sarda et al, 
2012). High-throughput sequencing technologies have rapidly 
increased our knowledge on DNA methylation in ecdysozoans 
(encompassing nematodes and arthropods), especially insects 
(Lyko et al, 2000; Elango et al, 2009; Park et al, 201 1; Zwier et al, 
2012), bringing new insights into the evolution of DNA methy- 
lation functions. Now, these questions have to be addressed in 
more distant invertebrate taxa such as lophotrochozoans (encom- 
passing worms and molluscs), which remain widely underde- 
scribed. Indeed, despite methylated DNA being present in the 
snail Biomphalaria glabrata (Fneich et al., 2013), the scallops 
Chlamys farreri and Patinopecten yessoensis (Wang et al., 2008), 
the truncated wedgeshell Donax trunculus (Petrovic et al, 2009) 
and the pacific oyster Crassostrea gigas (Gavery and Roberts, 2010; 
Riviere et al., 2013), only one single-base resolution methylome is 
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available to date in a mollusc, the pacific oyster Crassostrea gigas 
(Gavery and Roberts, 2013). 

C. gigas represents an interesting species for the study of DNA 
methylation in lophotrochozoans, with regards to its peculiar 
life traits, economic and ecological importance. Indeed, C. gigas, 
undergoes a pelagic development and metamorphoses before a 
benthic adult phase in the highly stressful intertidal area. Besides, 
oysters are successive hermaphrodites due to a yearly gonad 
renewal from stem cells. Therefore, the entire life cycle of oys- 
ters is punctuated with dramatic morpho-physiological changes, 
which rely on the implementation of transient transcriptomes 
within changing environments. The control of these transcrip- 
tomes likely implicates epigenetic mechanisms, which remain to 
be elucidated. The recent characterization of the C. gigas genome 
(Zhang et al., 2012) enabled an assessment that DNA methyla- 
tion was mostly intragenic in the oyster (Gavery and Roberts, 
2013). However, our recent results indicate that oyster methyla- 
tion patterns display temporal variations and could be uniquely 
biased toward the 5'-upstream region of gene subsets depending 
on physiological contexts (Riviere et al, 2013). Therefore, in this 
work we first review the literature on invertebrate DNA methyla- 
tion focusing on gene body methylation (GBM) and highlight the 
current knowledge in the oyster. Then, we discuss the perspec- 
tive of functional outcomes of possible specific DNA methylation 
features in C. gigas. 

DNA METHYLATION IN INVERTEBRATES 

EVOLUTION OF METHYLATION PAHERNS IN INVERTEBRATES 

The prevalence of DNA methylation within gene bodies in most 
invertebrates studied to date suggests that GBM is the ancestral 
state of this phenomenon (Sarda et al, 2012). Because methy- 
lated cytosines tend to spontaneously deaminate into thymines, 
an underrepresentation of the CpG dinucleotide in a region of 
interest reflects a sustained methylation (Bird, 1980). Therefore, 
the methylation level of genes is often inferred from the "nor- 
malized CpG content" (CpG O/E ratio), which compares the 
observed to the expected CpG dinucleotide contents (Shimizu 
et al, 1997). Several studies (Suzuki et al, 2007; Elango et al., 
2009; Gavery and Roberts, 2010) confirmed the robust negative 
correlation between the methylation level in the germline, i.e., 
possibly inherited across generations, and the CpG O/E ratio. 
Because of the increased mutation rate of methylated cytosines, 
DNA methylation is considered important in genome evolution. 
Accordingly, distant lineages present different patterns of GBM 
(Zemach et al., 2010; Sarda et al., 2012). Among divergent inver- 
tebrate groups, Sarda and colleagues reported that ecdysozoans 
(i.e., the honeybee Apis mellifera and the silkworm Bombyx mori) 
exhibit a lower GBM than the cnidarian anemone Nematostella 
vectensis and the tunicate sea squirt Ciona intestinalis (Sarda et al., 
2012). This indicates that protostomes (encompassing ecdyso- 
zoans and lophotrochozoans), have evolved toward a loss of (gene 
body) methylation, when compared to deuterostome animals 
(encompassing vertebrates) and to their common ancestor. The 
localization of methylated cytosines within transcription units 
also displays lineage-specific differences. Insects tend to have a 
higher methylation of the 5' than of the 3' regions of genes, 
whereas these levels are not different in the sea squirt, whUst their 



common ancestor, represented by the sea anemone, shows an 
intermediate profile (Zemach et al., 2010; Wang et al, 2013). The 
distinct evolutionary histories of GBM between taxa are further 
illustrated by the relationship between methylation and transpos- 
able elements (TE). Indeed, in contrast to vertebrates (Yoder et al., 
1997), methylation is clearly not the main mechanism used by 
insects in order to silence the genes lying within TEs (Lyko et al, 
2010; Wang et al., 2013), although ants demonstrate a high TE 
methylation (Bonasio et al, 2012). 

THE ROLE OF METHYLATION IN INVERTEBRATE GENOME EVOLUTION 

A general association between methylation of transcription units 
and protein conservation in invertebrates emerges from studies 
on insects (Hunt et al, 2010; Wang et al., 2013), cnidarian, and 
tunicates (Zemach et al, 2010; Flores et al., 2012; Sarda et al., 
2012). Indeed, despite slight particular features in the pea aphid 
(Hunt et al., 2010), highly methylated genes have a higher num- 
ber of orthologs than sparsely methylated genes (Sarda et al., 
2012). This is surprising because the hyper-mutability of methyl- 
cytosines would be expected to make DNA sequences diverge 
over evolutionary time (Cooper and Krawczak, 1989; Elango 
et al., 2008; Chuang and Chen, 2014). However, the weak cyto- 
sine representation within low CpG O/E genes could result in 
a lower frequency of methylation-dependent mutations (Hunt 
et al., 2010). Besides, mutations within genes methylated across 
many generations (low CpG O/E) are likely deleterious and thus 
would not be sustained in the genome. 

RELATIONSHIP BETWEEN DNA METHYLATION AND TRANSCRIPT 
VARIANTS IN INVERTEBRATES 

Despite distinct evolution histories between taxa, the role of GBM 
in transcript variant selection seems to be conserved. For exam- 
ple, GBM is widely implicated in the generation of alternative 
transcripts in mammals (Maunakea et al., 2010, 2013). Likewise, 
exon methylation is inversely correlated to exon skipping, and 
alternative splicing is increased among methylated genes in the 
honeybee (Flores et al., 2012). Consistently, Nasonia exons are 
"tagged" with cytosine methylation whereas introns are mostly 
unmethylated (Wang et al., 2013). Nevertheless, such patterns 
are not mandatory for a role of GBM in exon skipping. Indeed, 
Nasonia differentially spliced genes do not have an increased 
methylation probability, and methylated and non-methylated 
genes do not differ in their degree of alternative splicing (Wang 
et al., 2013). In the oyster, exons are preferentially methylated, 
albeit introns also exhibit significant methylation (Gavery and 
Roberts, 2013). Therefore, caution should be taken when infer- 
ring the role of methylation in exon selection from phyloge- 
netic proximity and/or methylation pattern similarity between 
species. 

GBM IN INVERTEBRATES IS ASSOCIATED WITH TRANSCRIPTIONAL 
REGULATION 

A greater number of transcripts can arise from methylated genes 
than from unmethylated genes in insects, possibly because methy- 
lation influences exon inclusion during transcription (Flores 
et al., 2012). In the oyster, a low inherited GBM (high CpG 
O/E) is supposed to increase the "transcriptional opportunities." 
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This refers to the association of DNA methylation with biases 
in the binding of the transcription proteins, thereby increas- 
ing the number of transcript variants (Roberts and Gavery, 

2012) . Thus, GBM would provide a mechanism for enhanced 
phenotypic plasticity from a limited number of genes. Parallel 
to transcript variant selection, the conservation of GBM also 
favors its biological significance in the regulation of transcrip- 
tion levels. In vertebrates, in contrast to promoter methyla- 
tion which generally represses transcription (Bird, 1992; Deaton 
et al, 2011; Vinson and Chatterjee, 2012), GBM is associated 
with active transcription (Ball et al, 2009; Maunakea et al., 
2010). In invertebrates which are mostly considered devoid of 
significant promoter CpG methylation (Elango and Yi, 2008), 
GBM also influences mRNA expression depending on the den- 
sity of methylation (Feng et al., 2010; Xiang et al., 2010; 
Zemach et al., 2010). Genome-wide observations across inver- 
tebrates indicate that transcriptional units are either highly 
or slightly methylated, and thereby display a general bimodal 
distribution regarding methylation. Nonetheless, the relation- 
ship between methylation and transcription exhibit species- 
specific features. Indeed, a positive linear correlation between 
GBM and mRNA levels exists in Nematostella vectensis and 
Bombyx mori (Xiang et al, 2010; Zemach et al., 2010). By 
contrast, in the honeybee most genes that are either highly 
or weakly expressed are heavily methylated (Zemach et al, 
2010). In adult oyster gUls, transcript levels increase with GBM 
until the 40th percentile and remain stable until the 100th 
percentile where methylation decreases (Gavery and Roberts, 

2013) . 

It is remarkable that gene subsets with distinct methylation 
levels correspond to distinct functions. In insects, genes with 
high canonical methylation (low CpG O/E) tend to be con- 
stitutively expressed at moderate to high levels and related to 
housekeeping functions, a conserved trend in the distant inver- 
tebrates investigated (Elango et al., 2009; Hunt et al., 2010; 
Sarda et al., 2012). Conversely, genes exhibiting little inher- 
ited methylation (high CpG O/E) display a broader range of 
mRNA levels and are implicated in regulated and/or inducible 
functions such as cell signaling and environmental stimulus 
response (Hunt et al, 2010). Accordingly, environmentally- 
induced methylomes control caste-related phenotypes in ants 
and honeybees (Lyko et al., 2010; Bonasio et al, 2012). In 
molluscs, similar interpretations were inferred from CpG O/E 
ratios (Gavery and Roberts, 2010; Fneich et al., 2013), further 
supported by the direct influence of DNA methylation on the 
transcription of development genes in the oyster (Riviere et al, 
2013). 

PERSPECTIVE: ARE THERE UNIQUE DNA METHYLATION 
FEATURES IN OYSTERS WITH FUNCTIONAL 
CONSEQUENCES? 

In molluscs, cytosines are more methylated than in other inverte- 
brate taxa (ca. 2 vs. ca. 0.15% in insects) (Zemach et al, 2010; 
Fneich et al, 2013; Gavery and Roberts, 2013). This indicates 
that molluscs exhibit specific characteristics in DNA methylation, 
suggesting more subtle roles and functional outcomes for this 
epigenetic mark than originally considered. 



IS THERE A "CGI PROMOTER-LIKE" TRANSCRIPTIONAL REGULATION 
IN OYSTERS? 

Intriguingly, a negative correlation exists between the mRNA level 
and the specific methylation of some oyster development genes 
(Riviere et al., 2013). This apparently contradicts the associa- 
tion of invertebrate GBM with active transcription (Ball et al., 
2009; Maunakea et al, 2010; Xiang et al, 2010; Zemach et al, 
2010; Gavery and Roberts, 2013), except if one considers that the 
methylated regions could extend upstream of gene bodies to the 
5' domain. Interestingly, such regions corresponding to the pre- 
dicted proximal promoter and/or first exon exhibit CpG islands 
(CGIs) (Figure 1). Due to the limitations in the MeDIP-qPCR 
methodology used, it remains unclear which of these regions 
are the actual targets of methylation. However, the recent single- 
base resolution methylome of adult C. gigas gills (Gavery and 
Roberts, 2013) reveals a significant methylation of putative pro- 
moter regions. Together, these observations suggest a bias in DNA 
methylation toward the 5' region of genes, which diminishes tran- 
scription at least for some genes at defined physiological states in 
the oyster. Therefore, we hypothesize that a "CGI promoter-like" 
transcriptional regulation could exist in C. gigas, in contrast with 
general findings in invertebrates (Elango and Yi, 2008; Zemach 
et al, 2010; Sarda et al., 2012), but in line with the results 
concerning the aplysia CREB2 gene (Rajasethupathy et al., 2012). 

DO DYNAMIC CHANGES IN DNA METHYLATION UNDERLY 
PHYSIOLOGICAL PLASTICITY? 

Consistent with variations in the amount of cytosine methylation 
during early oyster development (Riviere et al, 2013), our recent 
MeDIP-seq analyses indicate that developmental stages differ not 
only in the number of methylated genes, but also in the number 
of methylated regions per gene (Figure 2A). In adults, additional 
information comes from the study of the gametogenesis. Indeed, 
the methylation profiles are different between sexual resting (stage 
0) and maturity (stage 3) as well as between males and females 
(Figure 2B) (Riviere et al., in preparation). Triploid oysters were 
originally thought to be sterile because of aneuploidy-induced 
blocking of the gonial mitosis ("beta" triploids). However, some 
animals ("alpha" triploids) can escape sterility and produce 
mature gametes (Jouaux et al., 2010). Interestingly, methylation 
patterns in fertile alpha triploids are significantly different from 
those of sterile beta animals, yet are similar to those of fertile 
diploid oysters (Figure 2B). We found that specific methylation 
profiles are dynamic and match specific embryonic and germi- 
nal phenotypes in oysters. Hence, the association between DNA 
methylation and the regulation of differentiation genes in stem 
cells in mammals (review in Jones, 2012). Therefore, we assume 
a central role for DNA methylation in the differentiation of stem 
cells in C. gigas. 

Considering gene expression as a fundamentally random, 
stochastic phenomenon (Laforge et al, 2005; Raj and van 
Oudenaarden, 2008; Kupiec, 2010), DNA methylation would bias 
the number of probabilities in gene transcription in a given 
cellular context. As a consequence, methylomes would under- 
lie transcriptomes at the cellular level. A consistent hypoth- 
esis proposed GBM to increase the possibility of transcript 
variants from genes with low germline methylation, thereby 
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FIGURE 1 I The methylation of 5' regions can drive transcriptional 
repression during oyster development. Examples of development genes 
showing patterns suggestive of DNA methylation bias in upstream 
domains. DNA methylation (red) investigated by MeDIP-qPCR and cognate 
mRNA expression (blue) measured by RT-qPCR are shown for the C. gigas 
Abd, Alx4 and Post2 putative orthologs (development stages and GenBank 
accession numbers are indicated). The diagrams represent the genomic 
and methodological contexts. Scaffold numbers correspond to the 
fragment considered in the present assembly of the oyster genome (Zhang 
et al., 2012), start and stop positions (italic); intergenic (thin line), introns 
(thick line), exons (rectangles), putative MeDIPped fragments (green), qPCR 
amplified regions (orange), and CpG islands (red) are represented (modified 
from Riviere et al., 2013). *p < 0.05 Pearson or Spearman's correlation test 
between methylation and mRNA expression. 



underlying the phenotypic plasticity that allows C. gigas to adapt 
to unpredictable conditions at the population level (Roberts 
and Gavery, 2012). Therefore, the high CpG O/E ratio of oys- 
ter development genes, corresponding to an "inducible" profile, 
is inconsistent with their determined transcriptional patterns 
required in the context of a "fixed development program." In 
contrast, it rather fits with the present assumption that the devel- 
opment of C. gigas relies on the environmental stabilization of 
stochastic transcriptomes through epigenetic mechanisms like 
DNA methylation. Our hypothesis is supported by both the 



variation in the amount of methyl-DNA during embryogenesis 
and the suspected presence of promoter DNA methylation influ- 
encing oyster developmental gene transcription (Riviere et al., 
2013). From these observations which broaden former hypothe- 
ses (Roberts and Gavery, 2012; Gavery and Roberts, 2014), we 
propose that DNA methylation bears a crucial role in most if 
not all physiological processes which implicate the implemen- 
tation and the dynamic regulation of transient transcriptomes 
throughout the oyster life cycle. 

DOES DNA METHYLATION IN THE OYSTER BRIDGE THE GAP BETWEEN 
CELL DIFFERENTIATION, ENVIRONMENTAL ADAPTATION, AND 
EVOLUTIONARY SUCCESS? 

It remains to be elucidated whether methylation governs tran- 
scription, whether it just reflects local transcriptionaUy-active 
open chromatin, or both. Regardless, we assume that dynamic 
genome-wide DNA methylation profiles are integrated through 
the canalization of random interactions with the surrounding 
molecular machinery in the chromatin context (Badeaux and 
Shi, 2013) and lead to the implementation of environmentally- 
induced transcriptomes. Successful transcriptomes enable the 
development of sexually mature animals that may imprint methyl 
marks to DNA into their germinal cells, thereby passing suc- 
cessful methylomes to their progeny. If repeated across oyster 
generations, such successful methylomes could induce mutations 
and permanent changes in the genome. Therefore, in oysters 
DNA methylation could constitute a pivotal mechanism underly- 
ing processes from (i) cell differentiation during embryogenesis, 
(ii) organism survival in changing environments (Vandegehuchte 
and Janssen, 2013), and (iii) ultimately species adaptation over 
evolutionary time (Keller and Taylor, 2008). This confers a 
major role for DNA methylation in cell Darwinism and onto- 
phylogenesis in oysters (Maresca and Schwartz, 2006; RoUo, 
2006). However, one should remain cautious with this hypoth- 
esis because transgenerational imprinting has, to our knowledge, 
never been demonstrated in lophotrochozoans. 

FUTURE DIRECTIONS IN DNA METHYLATION RESEARCH IN 
OYSTERS 

The possible relevance of promoter methylation has been 
neglected in most invertebrates noticeably because of weak CpG 
representation and/or methylation in promoter regions (Elango 
and Yi, 2008). However, the transcriptional repression of the 
Aplysia CREB2 gene by the methylation of a CpG island in 
the promoter (Rajasethupathy et al., 2012) parallels our find- 
ing that 5' -methylation and mRNA expression can be negatively 
correlated (Riviere et al., 2013). Such a "CGI promoter-like" reg- 
ulation in invertebrates remains to be investigated, and in vitro 
methylated reporter gene constructs might be helpful. 

Current genome-wide studies bring informative static pic- 
tures of invertebrate methylomes. Forthcoming investigations 
should benefit from high-throughput technologies to decipher 
the dynamics of such methylomes and the mechanisms under- 
lying their functional implications. Future research should also 
address the complete epigenetic landscape of the chromatin 
beside cytosine methylation, including nucleosome occupancy 
and histone marks. Indeed, in flies (Ebert et al, 2006) and 
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FIGURE 2 I Dynamic changes in the exome-wide methylation 
pattern during the oyster development (A) and gametogenesis (B). 

Exon mapping after MeDIP-seq at the different stages indicated. 
Values are given as the mean ± s.e.m. of triplicate experiments. P 
values are given for One-Way ANOVA (p < 0.05 was considered 



mammals (Rach et al, 2011) similar DNA methylation levels 
have different outcomes, depending on the chromatin context. 
Furthermore, histone modifications are potentially relevant in 
oyster development and reproduction (Fellous et al., 2014). 

How the environment induces changes in DNA methyla- 
tion, especially in physiological contexts like development and 
reproduction, and whether these changes are persistent are also 
major issues. Indeed, as in mammals (Gabory et al, 2011), 
diet inputs during larval life control the expression of caste- 
related gene subsets in the honeybee (Kucharski et al., 2008). 
Moreover, investigations in Daphnia magna (Vandegehuchte 
et al., 2009, 2010) are suggestive of the trans-generational per- 
sistence of environmentally-induced epigenetic modifications in 
arthropods. Such "imprinting" mechanisms and the implication 
of the DNA methylation within remain to be demonstrated in 
lophotrochozoans, especially in "invasive" species like C. gigas. 
In this context, the estimated low inherited methylation of 
"inducible" oyster genes (Gavery and Roberts, 2010) suggests a 
possible reversibility of DNA methylation which remains to be 
explored in invertebrates. 

These questions open a major field of research and war- 
rant in-depth studies in lophotrochozoans, which encompass a 
tremendous variability of species and life traits. Many different 
reasons support the oyster as a model of interest in this perspec- 
tive: their peculiar ecology, the availability of their genome, the 
development of functional tools, the worldwide distribution of 
the species, its economic importance, and the present perspec- 
tive of functionally relevant specific DNA methylation features. 
A better knowledge of the dynamics, interactions, and roles of 
DNA methylation in Crassostrea gigas could greatly improve our 
understanding of this epigenetic mark in invertebrates. 
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